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A search for dinucleon decay into pions with the Super-Kamiokande detector has been performed 
with an exposure of 282.1 kiloton-years. Dinucleon decay is a process that violates baryon number 
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by two units. We present the first search for dinucleon decay to pions in a large water Cherenkov 
detector. The modes ^®0(pp) —>■ ^'’C7r"'"7r"'', ^®0(pn) —>■ and ^®0(nn) —>■ are 

investigated. No significant excess in the Super-Kamiokande data has been found, so a lower limit 
on the lifetime of the process per oxygen nucleus is determined. These limits are: > 

7.22 X 10®^ years, Tp„_>.„.+^o > 1.70 x 10^^ years, and T„„^^o.n.o > 4.04 x 10^^ years. The lower 
limits on each mode are about two orders of magnitude better than previous limits from searches 
for dinucleon decay in iron. 

PACS numbers: 13.30.Eg, 14.20.Dh, 29.40.Ka 


I. INTRODUCTION 

Baryon number seems to be a conserved quantity in or¬ 
der to explain the stability of matter, yet baryon number 
must be violated to explain the predominance of mat¬ 
ter over antimatter in the universe [T]. Single nucleon 
decay processes that violate baryon number by one unit 
but conserve A\B — L\ (where B is baryon number and 
L is lepton number) have been studied extensively at 
Super-Kamiokande w in the context of Grand Uni¬ 
fied Theories (GUTs) that predict an unstable proton. 
Lower limits on proton decay modes such as p —>■ e^7r° 
of 1.2 X 10^'* years [7] have ruled out the simplest GUTs, 
such as minimal SU (5) [5]. Other processes, such as those 
that violate B by two units, are of great theoretical inter¬ 
est although not necessarily predicted by GUT-inspired 
models. Searches for such processes constitute an im¬ 
portant subset of searches for baryon number violation, 
both in their own right and given the continued non¬ 
observation of single nucleon decay. Discovery of such a 
process would be a spectacular revelation of new physics, 
while non-observation provides important constraints on 
the theories which predict them. Dinucleon decay - the 
simultaneous decay of two bound nucleons into leptons 
or mesons - constitutes such a process. In this pa¬ 
per we search for the processes ^®0(pp) —>■ ^^G7r+7r+, 
—>■ ^^Ntt+tt^, and ^®0(nn) —)■ (abbre¬ 
viated throughout to pp —>■ pn —>■ and 

nn —>■ TT^TT^). 

Processes violating baryon number by two units can 
occur in a broad class of models [S]. The experimen¬ 
tally accessible processes of this kind are dinucleon decay 
and neutron-antineutron oscillation (n-nbar). The latter 
process may also be pursued in experiments with free 
neutrons uni, whereas dinucleon decay requires the close 
overlap of bound nucleons. The tree-level Feynman dia¬ 
gram for n-nbar and dinucleon decay to pions is the same, 
with two spectator quarks in the latter case. An exam¬ 
ple of a dinucleon decay diagram from one of the models 
discussed in [9] is shown in Fig. 1. The expected time 
scale for dinucleon decay to pions can be related to the 
time scale for n-nbar, if the same operator is dominant 
for each process. This is discussed further in Appendix 
A. 
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FIG. I. Feynman diagram for dinucleon decay to pions from 
one of the models discussed in [9], with spectator quarks 
gi ,?2 € {u,d). The Xi and X2 scalar particles belong to 
(6,1, —|) and (6,1, |) representations of the Standard Model 
SU{3) (g) SU{2) (g) U{1) gauge group. 

One dinucleon decay search has been performed be¬ 
fore at Super-Kamiokande, namely a search for decay 
to charged kaons m- The Frejus experiment, which 
used iron nuclei, performed an extensive search cover¬ 
ing dinucleon decay final states [T^]. The search in [T^] 
includes the modes analyzed in this paper. The more 
exclusive limits in this paper relative to those in |12) re¬ 
flect the large exposure of Super-Kamiokande. Our limit 
calculation accounts for systematic uncertainties, partic¬ 
ularly those due to intranuclear pion interactions. This is 
the first search for dinucleon decay into pions in a water 
Cherenkov detector. 

This paper proceeds as follows. Section II describes 
some basic features of the Super-Kamiokande exper¬ 
iment. Section III describes relevant aspects of our 
Monte-Carlo simulation, for both dinucleon decay sig¬ 
nal and atmospheric neutrino background. Section IV 
describes standard data selection and reconstruction al¬ 
gorithms used in Super-Kamiokande. Section V describes 
the analysis of all three modes, hrst discussing some fea¬ 
tures of boosted decision trees that are relevant for the 
pp —)■ 7r+7r+ and pn —>■ searches. Sections VI and 

VII describe our results and conclusion. 


II. THE SUPER-KAMIOKANDE DETECTOR 

Super-Kamiokande (hereafter “SK”) is a large, cylin¬ 
drical water Cherenkov detector. Full details of the de¬ 
tector can be found in m-: here, we note only the most 
basic features. 

The detector contains two optically-separated regions: 
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the inner detector (ID) and outer detector (OD). Inward¬ 
facing, 20-inch photomultiplier tubes (PMTs) collect 
Cherenkov light in the ID, while light from penetrating 
particles (typically cosmic-ray muons) is detected by 8- 
inch, outward-facing PMTs in the OD. The fiducial vol¬ 
ume of the ID is defined as the cylindrical volume with 
surfaces two meters inward from the ID wall. This cor¬ 
responds to 22.5 ktons, and 7.5 x 10^^ oxygen nuclei. 

SK has had four data-taking periods. It began taking 
data in April 1996, with 11,146 PMTs covering 40% of 
the ID surface. The first data-taking period, called SK-I, 
continued until July 2001, totaling 1489.2 livetime days, 
or 91.7 kton-years. An accident occurred in November 
2001, which caused the implosion of about half of the ID 
PMTs. The remaining 5,182 PMTs were re-distributed 
uniformly on the surface of the ID, providing 19% photo¬ 
coverage. To prevent such accidents in the future, each 
PMT has been enclosed in acrylic and fiber-reinforced 
plastic cases. The second data-taking period, which 
used this decreased photo-coverage, lasted from Decem¬ 
ber 2002 until October 2005, and is called SK-II. This 
lasted 798.6 livetime days, or 49.2 kton-years. Replace¬ 
ment PMTs were produced and the photo-coverage was 
returned to 40% in 2006. The third data-taking period, 
from July 2006 until September 2008, is called SK-III, 
comprising 518.1 livetime days (31.9 kton-years). The 
detector was upgraded in the summer of 2008, with im¬ 
proved electronics that record all PMT hit information 
without dead time m [IS]. This is the current con¬ 
figuration of the detector, called SK-IV. In this paper, 
SK-IV data are used up to March 2014, which comprises 
1775.6 livetime days (109.3 kton-years). The total SKI- 
IV dataset is 282.1 kton-years. 

III. SIMULATION 

Dinucleon decay signal and atmospheric neutrino back¬ 
ground Monte Carlo simulations (MC) are generated to 
estimate signal efficiency, expected background rates, 
and systematic uncertainties. Since each of the four 
SK data-taking periods have different detector configu¬ 
rations, separate sets of MC for both the dinucleon de¬ 
cay signal and the atmospheric neutrino background are 
used for each period. The dinucleon decay MC is based 
on single nucleon decay MC used for previous SK analy¬ 
ses HD. The atmospheric neutrino MC is the same MC 
sample used in other SK analyses, as described in detail 
in pre]. 

Dinucleon decay 

As dinucleon decay requires the close overlap of two 
nucleons, we restrict attention to the oxygen nucleus for 
the decaying nucleons. Fermi momentum, correlation be¬ 
tween nucleons, nuclear binding energy, and pion-nucleon 
interactions are taken into account. This last effect is es¬ 


pecially prominent, given the dense nuclear material the 
pions traverse. 

Fermi momentum and nuclear binding energy are 
taken from measured electron-^^C scattering [T7|. Nu¬ 
clear binding energy is taken into account by modifying 
the nucleon masses. Reference m estimates that 10% of 
single nucleons that undergo decay have wavefunctions 
that are correlated with other nucleons in the nucleus; 
we assume the same level of correlation for dinucleon de¬ 
cay pairs. Such correlated decay events typically have a 
lower total invariant mass, due to the momentum of the 
correlated nucleon, which is below Cherenkov threshold 
(Fig-i- The positions of the decaying nucleons in 
are approximated using the Woods-Saxon nuclear den¬ 
sity model m- Dinucleon decay leaves the remaining 
nucleus in an excited state from which it de-excites by 
the emission of gamma rays. However, the de-excitation 
gamma rays are on the order of a few MeV, while the pi¬ 
ons are on the order of a GeV, so the extra light coming 
from the gamma rays is negligible relative to the pions. 
The main difference in dinucleon decay simulation from 
single-nucleon decay is the available phase space. For din¬ 
ucleon decay, the initial total momentum and energy are 
summed over both nucleons (or three nucleons, if there 
is a correlated decay). 
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FIG. 2. Invariant mass of the pions from pp —>■ tt+tt"'", before 
nuclear effects. Plots for pn —> and nn —> 7r°7r° look 

similar. The dashed histogram shows the spectrum for events 
with a correlated nucleon. 

Pion-nucleon interactions in the nucleus are simulated 
using NEUT [20] • The interactions modeled include 
quasielastic scattering, in which the incoming and outgo¬ 
ing pion have the same charge, but there is some energy 
loss; charge exchange (e.g. 7r+ -|- n —7r° + p); absorp¬ 
tion (e.g. 7r+ -I- n —)■ p); and hadron production (e.g. 
7r+ -I- p —>■ 7r+ -I- 7r° -I- p). Quasielastic scattering is rel¬ 
evant at all pion momenta, as is charge exchange, al¬ 
though it has a lower overall cross-section. Absorption 
has a sizable cross-section in the pion momentum range 
100-300 MeV/c range, but is negligle elsewhere. Hadron 
production is relevant for pion momentum greater than 
500 MeV/c, and is therefore relevant for most pions im- 
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mediately after the decay. With these effects taken into 
account, less than half of the generated events in each 
mode escape the nucleus in the original pion configura¬ 
tion (Table |l]). 

After the pions leave the nucleus, their propagation 
and Cherenkov light emission is passed through a custom 
detector simulation based on GEANT3 [2T]. Charged pi¬ 
ons are susceptible to further secondary interactions in 
the water, while neutral pions decay before such interac¬ 
tions. Propagation of charged pions through the water 
below 500 MeV/c is simulated by custom code based on 
[22j ■ while charged pions above this momentum are han¬ 
dled by GCALOR P5] . 
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TABLE I. Number of pions for pp —>■ 7r’'"7r"'" (top), pn —>■ tt+tt® 
(middle), and nn —>■ (bottom) after nuclear effects. Per¬ 

centages indicate the percentage of events with the specihed 
number of tt^’s (columns) and number of tt^’s (rows). Entries 
with a “—” indicate a small contribution (<0.1%). The num¬ 
ber of pions prior to nuclear effects is shown in bold typeface. 


Monte carlo events are generated in a broadened fidu¬ 
cial volume (FV), greater than one meter from the ID 
wall as opposed to two meters which defines the FV. This 
allows us to account for events that migrate across the 
FV boundary when reconstructed. The final detection 
efficiency is defined as the ratio of events surviving all 
selection criteria to the total number of events generated 
in the fiducial volume. For each SK period, we generate a 
total of 25,000 events. Three times this number of events 
is generated in modes for which we use a multivariate 
search method. This is because the multivariate method 
requires three separate samples, as described later. 


Atmospheric neutrinos 

The atmospheric neutrino MC sample is based on the 
Honda atmospheric neutrino flux [24], and the NFUT 
neutrino interaction and nuclear effects simulator [20] . 
The dominant neutrino interaction background for the 
dinucleon decay searches comes from charged-current sin¬ 
gle pion production (CCItt, v + N N’ -\-tt ioi nucle¬ 

ons V, N' and charged lepton £), as well as both charged 
and neutral-current deep inelastic scattering (CCDIS and 
NCDIS, v+N —)■ l+N'+ hadrons, where i can be charged 
or neutral). Pions in CCItt interactions are produced 
following the model of Rein and Sehgal . This is the 
dominant hadron production interaction for hadronic in¬ 
variant mass W below about 2 GeVjf?. Deep inelastic 
scattering is treated differently depending on W. For 
1.3 GeVjf? < IT < 2.0 GeVjc^, only pions are consid¬ 
ered as outgoing hadrons, and the number of pions is 
two or greater to differentiate it from CCItt. For W > 
2.0 GeVjc^, other mesons are also considered, such as p 
mesons and kaons. After leaving the nucleus, particles 
are propagated through the water using the same tools 
described above for dinucleon decay. 

The atmospheric neutrino samples correspond to an 
exposure of 500 years for each SK period. Event rates 
are weighted to include the effect of neutrino oscillations. 
Final event rates after all selection criteria have been ap¬ 
plied are normalized to the corresponding SKI-IV detec¬ 
tor livetime. 


IV. REDUCTION AND RECONSTRUCTION 

This search uses the fully-contained (FC) dataset. 
Fully-contained events are defined as events in which 
all charged particles start and stop in the ID. Exiting 
and entering particles, particularly cosmic ray muons, 
are excluded by requiring minimal activity in the OD. 
The threshold for a triggered event in the FC dataset 
is 5.7 MeV of visible energy in SK-I, III, and IV, and 8 
MeV of visible energy in SK-II. Visible energy is defined 
as the energy to produce the observed light in the event 
if it were produced by a single electron. The majority 
of triggered events are cosmic ray muons and low energy 
radioactivity from material surrounding the detector. A 
dedicated set of event selection algorithms is applied to 
eliminate the majority of these backgrounds, such that 
all but a fraction of remaining events come from atmo¬ 
spheric neutrinos [16] . 

Events that remain are passed through a reconstruc¬ 
tion program. The reconstruction for each event de¬ 
termines an overall event vertex and the number of 
Cherenkov rings, as well as the direction, momentum, 
and particle classification of each ring. Events are as¬ 
sumed to originate from a single vertex, and the distribu¬ 
tion of observed charge is used to find the first Cherenkov 
ring. Additional rings are found using a Hough transform 
method [26] . Full details on our event reconstruction can 
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be found in m and P7] . 

Each Cherenkov ring is classified as “e-like” or “/i-like”, 
referring to electrons and muons. This is a conventional 
way of denoting whether a ring displays the characteris¬ 
tics of electromagnetic showers or not. Showering rings 
are produced by electrons and photons, but not protons, 
muons, or charged mesons. The classification of each 
ring is determined by its light pattern. Electromagnetic 
showers produce fuzzy rings due to scattering and pair 
production, while heavier particles produce rings with 
sharper edges. 

The momentum of each ring is determined based on the 
particle type, and the number of photo-electrons, referred 
to as charge, within a 70 degree half-angle cone of the ring 
direction. The charge is corrected for light attenuation 
and PMT angular acceptance. The overall charge in the 
cone is corrected for variations in the PMT gains. 

In addition to the standard ring-finding program, an 
algorithm described in detail in [55j is used to detect 
particle tracks. Unlike the standard ring-finder, this al¬ 
gorithm does not search for ring edges, but rather tests 
particle hypotheses as either electrons or photons (specif¬ 
ically, two photons from a tt^ decay) using only the light 
pattern. It has been used in atmospheric neutrino |29j 
and long-baseline [30] oscillation analyses to reject neu¬ 
tral current single tt® events. It is used in the pn —>■ 
search to identify a tt® from the decay and use its recon¬ 
structed mass as a discriminating variable for the analy¬ 
sis. 

After the primary event, electrons from stopping 
muons (Michel electrons) are found by searching for clus¬ 
ters of hits detected at the same time. The time window 
for such clusters extends to 20 ps after the main event. 
In the SK-I, II, and III data-taking periods, there was an 
impedance mismatch in the electronics which caused sig¬ 
nal reflection around 1000 ns after the main event. Thus, 
the time period 800-1200 ns after the main event was ex¬ 
cluded from the Michel electron search. The improved 
electronics of SK-IV eliminated this signal reflection, and 
as a result the tagging efhciency has improved from 80% 
to 96% for decays and 63% to 83% for decays be¬ 
tween SK-I-III and SK-IV. Hence, event selections that 
place an upper limit on the number of Michel electrons 
reduce the event rates more in SK-IV than in SK-I-III. 

After the reconstruction, a final set of selection criteria 
is applied to isolate fully-contained fiducial volume events 
(FCFV): 

• No significant OD activity: the largest OD hit clus¬ 
ter contains less than 10 hits (SK-I) or 16 hits (SK- 
II-IV). 

• The total visible energy is greater than 30 MeV. 

• The reconstructed event vertex must be at least 2 
m from the inner detector wall. This is the fiducial 
volume cut as defined in Section II. 

After the FCFV selection, the fraction of events coming 
from sources other than atmospheric neutrinos is esti¬ 


mated to be < 1%. In the next section, all analyses start 
with the FCFV events. 


V. ANALYSIS 

The experimental signature of a dinucleon decay to 
pions event, for all three modes considered, is two back- 
to-back pions, each forming at least one Cherenkov ring 
of the correct type. The characterization of such events 
is made difficult due to the interactions of charged pions 
in both the residual oxygen nucleus and the water. In¬ 
teractions in the water lead to truncated Cherenkov light 
signatures, which typically lead to poor momentum re¬ 
construction. Thus, for modes involving at least one tt'*' 
{pp —>■ and pn —>■ simple selection criteria 

based on kinematic variables such as the total momentum 
and invariant mass that are used in other nucleon decay 
analyses do not give optimal sensitivity. For these modes, 
a multivariate method is applied to statistically identify 
events with some, but perhaps not all, of the expected 
characteristics that differentiate signal and background. 
Specifically, a boosted decision tree (BDT) method is 
used. The method is implemented in the ROOT-based 
TMVA |3T] analysis library. This was also used in the 
dinucleon decay analysis in m- 

For nn —>■ the situation is more straightforward, 

as the 7r°’s that escape the nucleus decay to energetic 
photons before interaction in the water, and the electro¬ 
magnetic showers are well-reconstructed. For this mode, 
simple kinematic selection criteria are sufficient for the 
analysis. 


Boosted Decision Trees 

Here we describe some relevant details of the BDT, ap¬ 
plicable for both pp —>■ 7r''"7r+ and pn —> 7r"''7r'^ searches. 
Similar to other multivariate methods, a BDT takes a 
set of input variables - in this case, reconstructed infor¬ 
mation for both dinucleon decay signal and atmospheric 
neutrino background - with the goal of producing a single 
discriminating output variable that separates signal and 
background. To get to this single variable, three stages 
are required. They are referred to as training, testing, 
and analysis, and each uses separate sets of signal and 
background MC. We describe the process briefly. 

Boosted decision trees are sets of individual decision 
trees, simple binary trees formed by selection cuts on sin¬ 
gle variables. Trees are built using the training sample 
as follows. Starting with the full training sample, a selec¬ 
tion cut on a single variable is made that best separates 
signal and background. For the analyses in this paper, 
this optimal separation is determined by the Gini index 
g = p(l — p), where p is the purity of the sample after 
the selection. Other separation types exist, though they 
don’t lead to significantly different performance [3T| . The 
sample then splits off into two samples, or nodes, based 
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on this selection. Each of the daughter nodes does the 
same thing, separately finding the variable that best sep¬ 
arates between signal and background, and so forth until 
some stopping criterion is reached (such as a minimum 
number of events remaining on a node). Terminal nodes 
are declared signal-like or background-like based on the 
purity of signal events in the node. 

The training process continues over many such trees 
(typicially 1000 or so) using the same training sample, 
with the prescription that events in subsequent trees are 
re-weighted, or “boosted.” Several boosting algorithms 
exist, and are described in [3T]. Upon re-weighting, the 
best set of discriminating variables may change. Vari¬ 
ables can also be used more than once in a given tree. 

Once the training sample has been processed, addi¬ 
tional datasets may pass through the trained set of trees. 
On an event-by-event basis, the final BDT output is given 

by 

Ybdt = ^ Wihi (a;). (1) 

^^tree . ^ 

i—l 

Here Ntree is the size of the ensemble of trees, Wi is 
the weight of each tree in the ensemble, determined by 
the boosting method, hi is the “performance” of each 
event for the tree (whether it ended on a signal or 
background-like node) and x is the set of input variables. 
Variables are given a ranking, or “importance,” a nu¬ 
merical value that roughly corresponds to the number of 
times the variable is used in the training process relative 
to the other variables in the set. 

After the BDT is trained, two additional, statistically 
independent sets of signal and background MC are neces¬ 
sary. The test samples are used to determine the optimal 
boost method, as well as some additional parameters. 
The analysis samples are used to obtain the final sensitiv¬ 
ity. Unlike in the training and testing process, no infor¬ 
mation is given to the BDT regarding the composition of 
the analysis samples as either signal or background. Both 
the test and analysis samples are also important in de¬ 
tecting training bias in the BDT. This bias is a situation 
in which the BDT has been trained to model statistical 
fluctuations in the training sample, as opposed to mod¬ 
eling from general trends in the input variables. Such 
overtraining can lead to different performance between 
MC samples, even if they are statistically similar. 

Table [n] quantifies the signal and background MC sam¬ 
ples used for each stage. Separate trees are trained, 
tested, and analyzed for each SK period. Real data are 
processed through the trained BDTs as a final step. 

pp ->■ 7r+7r+ 

The interactions of 7r“'"’s in the nucleus and the water 
complicate the event signature for pp —>■ Tr+Tr^. A set of 
selection criteria is applied to keep pp —> tt+tt’'' events, 
and reduce the atmospheric neutrino background. The 



training 

testing 

analysis 

dinuc. 

decay 

25000 ev. 

25000 ev. 

25000 ev. 

atm-v 

~750,000 ev. 
(150 yr.) 

~750,000 ev. 
(150 yr.) 

~1,000,000 ev. 
(200 yr.) 


TABLE II. Independent MC sample sizes used for each stage 
of the BDT for both the pp —>■ and pn —>■ 7r'*’7r° searches, 

prior to any event selection. Dinucleon decay sample sizes 
show only the number of events, while the atmospheric neu¬ 
trino sample sizes show the approximate number of events 
and the equivalent number of years. The same statistics are 
used for SKI-IV. 


key characteristics are that the event has two “back-to- 
back” rings, as the ring-finding algorithm is still good at 
reconstructing the ring direction even with the poor mo¬ 
mentum resolution, and that the rings are identified as 
/i-like. Due to interactions in the water, the tt+’s can sim¬ 
ulate showering effects and therefore be mis-identified as 
e-like. However, the converse - in which rings produced 
by electrons or photons are misidentified as /r-like - is 
rare. The following selection criteria are applied: 

(Al) There is more than one Cherenkov ring. 

(A2) The two most energetic rings are /i-like. These are 
assumed to correspond to the 7r+ candidates. 

(A3) The angle between the two most energetic rings is 
greater than 120 degrees, a minimal requirement 
for “back-to-back” pions. 

(A4) The total visible energy is less than 1600 MeV, ap¬ 
proximately the maximum visible energy expected 
from a dinucleon decay event {2Mp — 2M^+), where 
Mp is the proton and M,r+ the pion mass. 



SK-I 

SK-II 

SK-III 

SK-IV 

Efr.(%) 

11.2±0.2 

10.5±0.2 

12.0T0.2 

12.1±0.2 

Bkg. 

33±0.9 

17T0.5 

13T0.4 

45±1.2 

data 

27 

14 

8 

43 


TABLE III. Efficiencies, expected backgrounds, and data 
events after selection criteria are applied in the pp —>■ tt+tt"'" 
search. Background rates are scaled to the appropriate SK 
period livetime. Statistical errors are shown. 


The signal efficiency, expected background, and data 
events remaining after the initial selection criteria is 
shown in Table HI After applying these criteria, we in¬ 
put the remaining events into the BDT. The following 
set of input variables is used, ordered by their relative 
importance, as described in the BDT section above: 


(al) The angle between the two most energetic /i-like 
rings. It is expected to peak around a perfect 
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SK-I 



120 130 140 150 160 170 180 

Angle (degrees) 


SK-II 



Angle (degrees) 


SK-III 



120 130 140 150 160 170 180 

Angle (degrees) 


SK-IV 



Angle (degrees) 


SK-I 



0 0.2 0.4 0.6 0.8 1 

Charge ratio 


SK-II 



SK-III 



0 0.2 0.4 0.6 0.8 1 

Charge ratio 



Charge ratio 


FIG. 3. The angle between ^-like rings and the ratio of charge carried by the most-energetic ring, for the pp —>■ 7r"''7r'^ search, for 
signal MC (solid histogram), atmospheric neutrino MC (dashed histogram), and data (crosses). Distributions are normalized 
to the number of events in the data. Details of each variable are described in the text. 






































































































































































“back-to-back” angle of 180 degrees for the signal, 
and have a falling spectrum for increasing angle in 
the background. 

(a2) The ratio of charge carried by the highest-energy 
ring to the total charge of all rings. This is ex¬ 
pected to be smaller for signal, since the pions are 
expected to have roughly equal energy. By con¬ 
trast, we expect most of the energy from a charged- 
current (CC) atmospheric neutrino event to be in 
the charged lepton ring. 

(a3) The total visible energy. This has a broad distri¬ 
bution up to 1600 MeV for the signal, but peaks at 
a lower value for the background around 200 MeV. 

(a4) The maximum distance between a Michel electron 
vertex and the primary event vertex. We expect 
this to be larger for energetic muons than for pi¬ 
ons, since the muons penetrate water without sud¬ 
den hadronic interactions. Muons above Cherenkov 
threshold sometimes appear in the signal (about 
9% of the time), but are much more common in 
the background, due to charged-current interac¬ 
tions such as single pion production. 

(a5) The maximum angle between any /r-like ring direc¬ 
tion vector, and any Michel electron vertex vector 
(as pointed to from the primary vertex). Similar 
to (al), this is designed to look for “back-to-back” 
events, and is an additional check on the topology 
of the /i-like rings. 

(a6) The magnitude of the vector sum of corrected 
charge associated with all rings, where the correc¬ 
tions are applied for various factors, as described in 
Section IV. This is similar to the total momentum, 
since the corrected charge associated with each ring 
is used to determine its momentum. Due to typi¬ 
cally poor momentum resolution from hadronic en¬ 
ergy losses, it is preferable to use this variable in¬ 
stead of the total momentum. 

(a7) The number of Cherenkov rings. Though we are 
trying to isolate two pions, they typically interact in 
the water, often producing visible secondary tracks, 
so that the number of rings distribution peaks at 
three. By contrast, the Cherenkov ring spectrum 
for the background falls with increasing ring num¬ 
ber. 

(a8) The number of Michel electrons. Naively we would 
expect two Michel electrons from two charged pi¬ 
ons. However, the pions could undergo charge ex¬ 
change or absorption, leading to zero or one Michel 
electrons. Though the signal and background dis¬ 
tributions both peak at one decay electron, this 
variable is relatively independent of most of the 
other variables, which depend on Cherenkov ring¬ 
counting. 

(a9) The number of ^-like rings. Similar to the to¬ 
tal number of rings, this variable specifically seeks 
scatters that produce a secondary charged pion. 


The relative importance of all input variables is shown 
in Table IV Fig. shows plots of the two most discrim¬ 


inating variables, the angle between tt’’' candidates and 
the ratio of charge carried by the most energetic ring. 
The number of events in the data is small after the se¬ 
lection criteria is applied, but the data and atmospheric 
neutrino MC are consistent. 

The final distributions of the BDT output for all SK 
periods are shown in Fig. For pp —> tt+tt''', the gra¬ 
dient boost [3T] method is used. The final efhciendes, 
backgrounds, and candidate data events are shown in 
Table |V] SK-II has a slightly lower efficiency than the 
other SK periods, as fewer events pass selection criterion 
(A2). 


Variable 

Importance 

Angle between /r-like rings 

0.16 

Ratio of charge carried by 

0.15 

most energetic ring 

Visible energy 

0.15 

Max. distance to 

Michel vertex 

0.13 

Max. angle between /r-like 
ring and Michel vertex 

0.13 

Magnitude of vector sum 
of corrected charge 

0.12 

Number of rings 

0.071 

Number of Michel electrons 

0.055 

Number of /r-like rings 

0.045 


TABLE IV. Relative importance of each variable in the 
pp —> search, averaged across SK periods. 



SK-I 

SK-II 

SK-III 

SK-IV 

Eff. (%) 

6.1T0.2 

5.3±0.2 

6.4T0.2 

5.8±0.2 

Bkg. (MT-yr) 

17.8±1.8 

14.3±1.6 

17.4±1.7 

14.2T1.6 

Bkg. (SK live.) 

1.6 

0.70 

0.56 

1.6 

Candidates 

0 

1 

0 

1 


TABLE V. Efficiency, expected background events, and can¬ 
didate data events for the SKI-IV pp —>■ tt+tt"'" search. Statis¬ 
tical errors are quoted for efficiency and background. Back¬ 
ground is quoted both for the appropriate SK livetime, and 
per megaton-year. 


mode 

SK-I 

SK-II 

SK-III 

SK-IV 

CCItt 

45±7% 

42±7% 

51±7% 

42±7% 

CCDIS 

30±5% 

31±6% 

24±5% 

34±7% 

NCDIS 

10±4% 

9±4% 

6±3% 

13±5% 

CCQE 

9±3% 

9±3% 

10±3% 

10±4% 

NCItt 

5±3% 

8±4% 

8±3% 

0+1% 


TABLE VI. Atmospheric neutrino interaction modes in the 
remaining background for the pp —> Tr^Tr"*" search by percent¬ 
age, as estimated by atmospheric neutrino MC. 


The final cut on the BDT output was determined using 
only signal and background MC, not the real data. The 
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SK-I 



BDT Output 



BDT Output 


SK-I I 



BDT Output 


SK-IV 



BDT Output 


FIG. 4. BDT output for pp —>■ 7r"'"7r^ signal MC (solid), atmospheric neutrino background MC (dashed), and data (crosses). 
The vertical green line indicates the BDT cut value, and the arrow indicates that only events to the right of the cut are kept. 


Variable 

Candidate 
SK-II SK-IV 

Signal 

SK-II SK-IV 

Background 
SK-II SK-IV 

Angle between /r-like rings (degrees) 

177.7 

151.8 

152.3 

152.1 

138.3 

138.1 

Ratio of charge carried by most energetic ring (dimensionless) 

0.60 

0.46 

0.55 

0.55 

0.66 

0.66 

Visible energy (MeV) 

400.0 

745.6 

680.5 

693.7 

440.2 

407.3 

Max. distance to Michel vertex (cm) 

1245.6 

1420.1 

440.7 

424.8 

523.0 

459.2 

Max. angle between /i-like ring and Michel vertex (degrees) 

166.3 

111.9 

94.2 

109.0 

93.2 

110.9 

Magnitude of vector sum of corrected charge (corrected p.e.) 

306.9 

1799.5 

601.3 

1310.0 

448.7 

797.5 

Number of rings (dimensionless) 

3 

3 

2.89 

2.95 

2.39 

2.40 

Number of Michel electrons (dimensionless) 

1 

2 

0.91 

1.09 

1.00 

1.21 

Number of /r-like rings (dimensionless) 

2 

2 

2.14 

2.16 

2.05 

2.05 


TABLE VII. Candidate variable values, and mean variable values for signal and background distributions, for each of the 
variables used in the BDT for the pp —>■ tt+tt''' search. Values are given for SK-II and SK-IV, for which single candidate events 


were found. 


cut takes two basic considerations into account: (I) the 
cut maximizes elVb, where e is signal efficiency and b is 
the background rate; and (2) the cut is not too close to 
the signal peak, since variations are propagated through 
the BDT to estimate systematics, and the part of the 
distribution near the peak changes the most under these 
variations. For (2), 0.7 was the closest value to the peak 
that we considered. As the efficiency (background) dis¬ 
tribution increases (decreases) nearly monotonically for 


all SK periods, it was found that a cut at 0.7 satisfied (1) 
while obeying the constraint of (2). 

The largest component (^40%) of the remaining at¬ 
mospheric background is from charged current single 
pion production (CCItt). In this mode, typically a 

will produce a muon and a charged pion at a 
large angle, thus passing the candidate angle re¬ 
quirement. Other contributing background interac¬ 
tion modes include charged-current deep inelastic scat- 
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FIG. 5. A pp —>■ 7r"'"7r^ MC event. The intensity scale indi¬ 
cates the amount of charge deposited on the PMTs. The main 
display shows the ID, while the small display on the top right 
shows the OD. The white -l-’s indicate the horizontal and ver¬ 
tical location of the reconstructed vertex, projected onto the 
detector wall. Orange rings are fit as p-like, and cyan rings as 
e-like. Solid rings indicate that the ring’s opening angle is near 
the ultra-relativistic limit in water of 6^ = cos~^ {^/n) « 41 
degrees, while dashed or dotted rings are more collapsed. The 
two solid orange rings correspond to true tt^’s, and have an 
angle of 164 degrees between them. The dotted cyan ring is 
fit as e-like, but is actually a hard scatter of one of the tt^’s. 



FIG. 6. Surviving atmospheric nentrino MG event for the 
pp —> search. The interaction type is charged- 

current single TT^ production. The uppermost p-like ring cor¬ 
responds to a p“, while the smaller p-like ring is a tt"*". The 
e-like ring is mis-fit, as it comes from a hard scatter of the 
tt"*". The angle between the two p-like rings is 170 degrees, 
thereby mimicking the main event signature. The momentum 
of the muon p-like ring is 402 MeV/c, and the momentum of 
the tt"*" p-like ring is 301.3 MeV/c 


tering (CCDIS), neutral-current deep inelastic scatter¬ 
ing (NCDIS), charged-current quasi-elastic (CCQE) and 



FIG. 7. pp —> TT^Tr"*" candidate event in SK-II. The two p- 
like rings are strikingly back-to-back, at almost exactly 180 
degrees to each other. The e-like ring could be due to a hard 
scatter. There is one Michel electron, not shown. The p-like 
rings have momenta of 399.0 MeV/c (solid) and 280.5 MeV/c 
(dashed), and the e-like ring has momentum 42.0 MeV/c. 
The visible energy is 387.0 MeV, the total momentum 124.2 
MeV/c, and the invariant mass 710.9 MeV/c^. 



FIG. 8. pp —>■ 7r’'"7r"'" candidate event in SK-IV. The two p- 
like rings are at an angle of 151 degrees to each other, with 
momenta of 541.4 and 389.5 MeV/c. The dashed e-like ring 
could be due to a hard scatter. There are 2 Michel electrons, 
not shown. The visible energy is 742.8 MeV, the total mo¬ 
mentum is 445.8 MeV/c, and the invariant mass is 1022.3 
MeV/c^ 


neutral-current single pion production (NCItt). Table VI 


quantifies the neutrino interaction modes of the remain¬ 
ing background. The percentages agree within statistical 
fluctuations across SK periods. 

A pp ^ 7r+7r+ MC event is shown in Fig. An at¬ 
mospheric neutrino event surviving all selection criteria. 
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including the final BDT cut, is shown in Fig. [6| Two 
candidate events were found, one in SK-II (Fig. and 
another in SK-IV (Fig.[^. No events are observed in SK-I 
or SK-III. Table [Vn| shows values of the BDT input vari¬ 
ables for each candidate, compared to the means of the 
signal and background distributions. While both candi¬ 
date events appear signal-like in some variables, partic¬ 
ularly the angle between /r-like rings (both candidates) 
and energy variables (SK-IV candidate), the candidates 
are consistent with the expected background rate for each 
SK period (Table [V|, as well as the total background rate 
of 4.5 expected across SKI-IV. Thus, we conclude that 
there is no evidence for pp —>■ tt+tt"'' in the data. 

pn 7r“*“7r'’ 

There are many similarities between the pn —>■ 
search and thspp —>■ tt+tt’'" search. The main difference is 
that the tt® provides a new set of discriminatory variables, 
and some of the multivariate input variables only apply 
to the reconstructed 7r°. 

First, we apply a set of selection criteria: 


Table VIII About a third of the primary 7r*^’s in the 
pn —>■ TT+TT® events surviving the selection criteria are 
reconstructed as a single e-like ring. This is mainly be¬ 
cause the typically high momentum of the 7r°’s leads to 
a small opening angle between the photons in the subse¬ 
quent 7r° —77 decay, so that the rings are not distinct. 

Events surviving the selection criteria are passed into a 
BDT. The following set of input variables is used, ordered 
according to variable importance: 


(bl) The momentum of the 7r° candidate. The 7r° mo¬ 
mentum peaks near 1 GeV for the signal, which 
is typically well-reconstructed. The background 
peaks at a much lower value of ^ 200 MeV. 

(b2) The angle between the reconstructed 7r+ and 7r°. 
This peaks near 180 degrees for the signal. 

(b3) The momentum of the candidate. The momen¬ 
tum resolution for charged pions from dinucleon de¬ 
cay is typically poor. Despite this, it still peaks at a 
higher value in the signal than the background (470 
MeV/c vs. 390 MeV/c), and is relatively indepen¬ 
dent of the better-reconstructed momentum. 


(Bl) There is more than one Cherenkov ring. 

(B2) There is at least one e-like and one /r-like ring. The 
7r+ is assumed to correspond to the most energetic 
/i-like ring, while the 7r° is assumed to correspond 
to the e-like ring (if one e-like ring is found) or the 
two most energetic e-like rings (if more than one is 
found). 

(B3) The number of Michel electrons is no more than 
one. 


(b4) The invariant mass of the 7r°. This can be esti¬ 
mated for both one and two-ring hypotheses, by 
means of the specialized algorithm referred to in 
Section ITVl 

(b5) The ratio of charge carried by the highest-energy 
ring to the total charge of all rings. 

(b6) The total visible energy. 

(b7) The number of Michel electrons. 


(B4) The “reduced” visible energy, defined as the total 
visible energy minus the energy of the as defined 

in the previous criterion, is less than 800 MeV. This 
is similar to criterion (A4) used in the pp —> tt+tt"'' 
search - the maximum visible energy from the tt"*' 
is Mproton - M^+ « 800 MeV. 

(B5) The angle between the reconstructed 7r+ and 7r° is 
greater than 120 degrees. 



SK-I 

SK-II 

SK-III 

SK-IV 

Eff. (%) 

21.0±0.3 

21.9±0.3 

21.6±0.3 

21.1±0.3 

Bkg. 

132±1.8 

69±1.0 

48±0.6 

147±2.0 

data 

136 

66 

45 

171 


TABLE VIII. Efficiencies, expected backgrounds, and data 
events after selection criteria in the pn —>■ search. Back¬ 

ground rates are scaled to the appropriate SK period livetime. 
Statistical errors are shown. 


The signal efficiency, expected background, and data 
events remaining after the selection criteria is shown in 


Variables (b5), (b6), and (b7) were described for the 
pp —>■ 7r+7r+ search. The ranking of these variables is 
shown in Table |IX| The two most discriminating input 
variables are shown in Fig. The 7r°-7r+ angle is highly 
discriminating for the same reason as the angle between 
7r+ candidates in the pp 7r+7r+ search. The 7r° mo¬ 
mentum is also a highly discriminating variable, since 
the peak near 1 GeV/c for the signal tt® momentum, ex¬ 
pected for back-to-back pions emitted with initial energy 
~ 2 GeV, is absent in the background. The data and 
atmospheric neutrino MG are in good agreement. 

The final BDT outputs are shown in Fig. 10 For the 
pn —>■ 7r+7r° search, an adaptive boost method [3T] was 
chosen. 

The final cut on the BDT output was designed so that 
less than one background event was expected across all 
SK periods, with contributions from each SK period or¬ 
dered by their livetime. It is determined with MG only. 
The total expected background is 0.75 events, consis¬ 
tent with the one candidate event found in SK-I. Unlike 
the pp —>■ 7r“'"7r^ search, the signal distributions are not 
sharply peaked, so that the efficiency does not change 
drastically as the cut value is shifted. 
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FIG. 9. The angle and 7r° momentum for the pn —>■ 7r"''7r° search. Signal and atmpospheric neutrino MC (solid and 

dashed histograms) are normalized to the number of events in the data (crosses) for each SK period. Details of each variable 
are described in the text. 
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SK-I 



BDT Output 

SK-III 



- 0.6 - 0.4 - 0.2 0 0.2 0.4 0.6 

BDT Output 


SK-I I 



BDT Output 


SK-IV 



- 0.6 - 0.4 - 0.2 0 0.2 0.4 0.6 


BDT Output 


FIG. 10. BDT output for pn —> 7r"'"7r° signal MC (solid), atmospheric neutrino background MC (dashed), and data (crosses). 
The green vertical line indicates the BDT cut value, and the arrow indicates that only events to the right of the cut are kept. 
One candidate event can be seen in the SK-I data distribution, just to the right of the cut at 0.19. 


Variable 

Importance 

TT*^ candidate momentum 

0.19 

Angle between 7r° and 

7r+ candidates 

0.17 

tt"*' candidate momentum 

0.16 

7r° candidate invariant mass 

0.15 

Ratio of charge carried by 

0.14 

most energetic ring 

Visible energy 

0.14 

Number of Michel electrons 

0.058 


TABLE IX. Relative importance of each variable in the 
pn —>■ Tv'^'K^ search, averaged across SK periods. 



SK-I 

SK-II 

SK-III 

SK-IV 

Cut 

0.19 

0.24 

0.20 

0.17 

Eff. (%) 

10.2±0.2 

10.0T0.2 

9.4T0.2 

10.4T 0.2 

Bkg. (MT-yr) 

2.7T0.7 

2.3T0.7 

2.2T0.7 

2.9T0.8 

Bkg. (SK live.) 

0.25 

0.11 

0.07 

0.32 

Candidates 

1 

0 

0 

0 


TABLE X. Final BDT cut value, efficiency, expected back¬ 
ground, and candidate data events for the SKI-IV pn —>■ tt+tt® 
search, with statistical uncertainties shown for efficiency and 
expected background. Background rates are quoted for the 
appropriate SK livetime, and per megaton-year. 


The cut on the final BDT output, efficiencies, back¬ 
ground rates, and candidates are shown in Table 
CCDIS and CCItt are the main components of the re¬ 
maining background, detailed in Table |XI| 

An example pn —>■ 7r“'"7r'^ MC event display is shown in 
Fig.[ni and an example background MC event is shown 
in Fig. ^ One candidate is found in SK-I, shown in 
Fig. Table |XII| shows the values of the BDT in¬ 
put variables for the candidate event, compared with 


the means of the signal and background distributions. 
While the candidate event appears closer to the signal 
mean values for energy-related variables, it appears more 
background-like for the 7r° — 7r+ angle and the recon¬ 
structed 7r° mass. There are no candidate events in SKII- 
IV data. Although the expected background in SK-I is 
much lower than one event, it is consistent within about 
Icr with the expected background of 0.75 events in the 
entire SKI-IV dataset. The probability of observing one 
or more background event in SK-I, without incorporating 
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FIG. 11. A pn —> 7r"'"7r° MC event. The /i-like ring corre¬ 
sponds to a true . The e-like ring contains two overlapping 
7 ’s from a 7r°. The angle between the p- and e-like rings is 
165 degrees. The e-like ring has a reconstructed momentum of 
893 MeV/c, while the original 7r° had momentum 876 MeV/c. 



FIG. 12. A surviving atmospheric neutrino MC event for 
the pn —>■ search. The interaction type is CCItt®, 

—>■ The two overlapping e-like rings are pho¬ 

tons from a 7r°, with total momentum 939 MeV/c. The /i-like 
ring corresponds to the , with momentum 410 MeV/c. The 
angle between them is 151 degrees, thereby appearing back- 
to-back. 


systematic uncertainties, is 22%, and therefore the event 
is consistent with background expectations. Further, the 
low 7r° mass of the candidate event, as reconstructed by 
the specialized algorithm discussed in Section |IVj sug¬ 
gests that it is not a pn ^ 7r+7r° signal. The e-like ring 
is more likely to be an electron from, for instance, a CCt'e 
Itt event. Thus, we conclude that there is no evidence 
for pn —)• 7r+7r° in the SK-I-IV dataset. 



FIG. 13. The pn —>■ candidate in SK-I. The e-like and 

/i-like rings are at an angle of 140 degrees. The e-like ring 
has a high momentum of 987 MeV/c, but its fit 7r° mass is 
low, at 10 MeV/c^. The /t-like ring has a momentum of 460.0 
MeV/c. The total momentum and invariant mass are 700.1 
MeV/c and 1266.4 MeV/c^, respectively. 


mode 

SK-I 

SK-II 

SK-III 

SK-IV 

CGItt 

33±17% 

29±17% 

43±19% 

45±18% 

GGDIS 

38±17% 

46±19% 

31±18% 

40±16% 

NGDIS 

8±8% 

20±14% 

21±15% 

8±8% 

GGQE 

8±8% 

5±5% 

5±5% 

8±8% 


TABLE XL Neutrino interaction mode percentages for re¬ 
maining background in the pn —>■ search for SKI-IV. 

The percentages agree across SK periods within statistical 
fluctuations, which are large due to the low number of re¬ 
maining events. 


Variable 

Candidate Signal Background 

7r“ momentum (MeV/c) 

986.5 

756.0 

233.7 

7r°-7r^ angle (degrees) 

142.8 

152.9 

140.4 

tt'*’ momentum (MeV/c) 

459.5 

470.3 

394.5 

7r° mass (MeV/c^) 

10.0 

143.0 

102.2 

Charge ratio (dimensionless) 

0.77 

0.56 

0.58 

Visible energy (MeV) 

1171.6 

1098.6 

495.0 

Michel electrons (dimensionless) 

0 

0.40 

0.51 


TABLE XII. Gandidate variable values, and mean variable 
values for signal and background distributions, for each of 
the variables used in the BDT for the pn —>■ tt+tt® search. 
Letters for each variable are assigned in the text. Values are 
given for SK-I, for which a single candidate event was found. 


The nn search is distinct from the other two 

searches, in that no multivariate method is used. Similar 
to other nucleon decay searches at SK Hi], the total 
momentum and invariant mass are highly discriminatory 
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variables, so that the search proceeds with a simple set of 
selection criteria. The criteria are similar to those used 
for the p —>■ e+7r° analysis [1]. 

(Cl) 2 < number of Cherenkov rings < 4; 

(C2) all rings are e-like; 

(C3) no Michel electrons; 

(C4) Ptot < 600 MeV/c; 

(C5) imQMeV/c^ < < 2000 MeC/c^. 



SK-I 

SK-II 

SK-III 

SK-IV 

FCFV 





Eff. (%) 

94.1T0.7 

96.6±0.7 

94.4±0.7 

94.1T0.7 

Bkg. 

11846T11.1 

6389±5.8 

4187±3.7 

14214T12.7 

Data 

12299 

6610 

4355 

14444 

(Cl) 





Eff. (%) 

76.7T0.6 

76.4T0.6 

76.6T0.6 

76.8T0.6 

Bkg. 

3434±5.6 

1846±3.0 

1204T2.0 

4135±6.7 

Data 

3558 

1969 

1239 

4215 

(C2) 





Eff. (%) 

61.0T0.5 

59.1T0.5 

59.6T0.5 

60.1T0.5 

Bkg. 

2131T4.3 

1175±2.3 

728±1.5 

2502±5.0 

Data 

2249 

1268 

782 

2623 

(C3) 





Eff. (%) 

54.4±0.5 

52.6T0.5 

52.8±0.5 

52.1T0.5 

Bkg. 

1441±3.4 

799±1.9 

497±1.2 

1543±1.1 

Data 

1492 

837 

529 

1627 

(C4) 





Eff. (%) 

30.1T0.4 

27.0T0.4 

29.1T0.4 

29.3T0.4 

Bkg. 

837±2.6 

470T1.4 

290±0.9 

940±3.0 

Data 

853 

503 

310 

1021 

(C5) 





Eff. (%) 

22.1T0.3 

18.8T0.3 

20.9T 0.3 

21.4T 0.3 

Bkg. 

0.05±0.02 

0.04T0.01 

0.03±0.01 

0.02±0.01 

Data 

0 

0 

0 

0 


TABLE XIII. Efficiency, expected background, and data 
events for the SKI-IV nn —>■ search, for each analysis 

selection criterion. Statistical errors are shown. The criteria 
are described in the text. 

The total momentum Ptot is defined as Ptot = 
where firing is the number of reconstructed 
rings, Pi is the 3-momentum of the ring, and the in¬ 
variant mass is defined as Minv = where 

Etot = Ip) I (the rings are assumed to come from 

photons, so there is no rest-mass energy). A nn — f 
MC event is shown in Fig. and a background MC 
event surviving all selection criteria is shown in Fig. |15| 
Criterion (Cl) allows for both one- and two-ring 7r° fits. 
Criteria (C2) and (C3) are designed to eliminate muons 



FIG. 14. A nn —>■ MC event. The two overlapping e-like 

rings on the top left correspond to two y’s from a tt® decay, 
and the other ring also corresponds to a 7r°. The trne mo¬ 
mentum of the 7r° corresponding to the two-ring fit is 872.7 
MeV/c, and the true momentum of the tt® corresponding to 
the one-ring fit is 936.7 MeV/c. The total reconstructed mo¬ 
mentum of this event is 219 MeV/c, and the reconstructed 
invariant mass is 1793 MeV/c^. 



FIG. 15. Surviving atmospheric neutrino MC event for the 
nn —>■ 7r°7r° search. The interaction is neutral current deep- 
inelastic scattering. In particular, the interaction produced a 
7T° and an rj meson, from which four y’s are produced. All but 
the left-most ring correspond to true photons from the 7r° and 
T] meson. Many additional charged pions are also produced, 
concentrated near the center of the display, thereby compli¬ 
cating reconstruction. The left-most ring corresponds to a 
charged pion. There is a ring near the center from multiple, 
closely overlapping charged pion tracks, which is not found 
by the reconstruction program. The total reconstructed mo¬ 
mentum is 553 MeV/c, and the reconstructed invariant mass 
is 1763 MeV/c^ 









































16 


SK-I 
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SK-I I 



Total Momentum (MeV/c) 



0 500 1000 1500 2000 2500 


Invariant Mass (MeV/c^) 


SK-III 



Total Momentum (MeV/c) 



SK-IV 



SK-IV 



0 500 1000 1500 2000 2500 
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FIG. 16. Total momentum (left) and invariant mass (right) for data (crosses), nn —>■ MC (solid histogram), and at¬ 

mospheric neutrino MC (dashed histogram). Criteria (C1)-(C3) have been applied to all distributions. The invariant mass 
criterion (C5) has been applied to the total momentum, and the total momentum criterion (C4) to the invariant mass. The 
thin green vertical lines indicate the locations of the final cuts for each variable. 
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FIG. 17. The total momentum and invariant mass distributions for nn —>■ MC (left), atmospheric neutrino MC (center), 

and data (right) for SK 1-IV. As seen in Fig. 16 there is no major structural difference for SK-I-IV distributions, so events from 
all periods are combined. Events have passed criteria (C1)-(C3). The solid box indicates the total momentum and invariant 
mass cuts ((C4) and (C5) in the text). 


and charged pions. Criterion (C4) is chosen to ensure a 
low total momentum, and to reject cases where only one 
TT^ is reconstructed, as each from nn —> has a 

momentum that peaks closer to 1 GeV/c. Criterion (C5) 
is chosen to cover the ^ 2M„ region as tightly as possible 
without getting too close to the peak (Fig. [^. This 
eliminates nearly all background, which peaks sharply 
near M^o. 


The signal efficiency, background, and data events af¬ 
ter each selection criteria are shown in Table IXIIII The 
final efficiencies and background rates are summarized in 
Table [XIV| The remaining background is mostly neutral 
current deep-inelastic scattering (Table XV). The inter¬ 
action mode percentages agree across SK periods within 
statistical fluctuations, which are large due to the very 
low number of remaining events. The lower background 
in SK-IV is due to improved Michel electron tagging rel¬ 
ative to SK-I, II, and III. 



SK-I 

SK-II 

SK-III 

SK-IV 

Efficiency (%) 

22.1±0.3 

18.8±0.3 

20.9± 0.3 

21.4± 0.3 

bkg. (MT-yr.) 

0.6±0.2 

0.9±0.3 

1.0±0.3 

0.2±0.1 

bkg. (SK live.) 

0.05 

0.04 

0.03 

0.02 

candidates 

0 

0 

0 

0 


TABLE XIV. Efficiency, expected background, and candidate 
data events for the SKI-IV nn —>■ TT^TT^ search, with statisti¬ 
cal uncertainties. Background is quoted both for each SK 
livetime and in megaton-years. 


Fig. shows the final total momentum and invariant 
mass distributions for signal and background MC, and 
the data. No candidate events are found for this mode, 
consistent with the near-zero expected background for all 
SK periods. 


mode 

SK-I 

SK-II 

SK-III 

SK-IV 

NCDIS 

63±32% 

30±18% 

67±25% 

50±50% 

CCDIS 

15±16% 

50±22% 

24±14% 

0+50% 

CCItt 

21±15% 

20±14% 

9±9% 

51+51% 


TABLE XV. neutrino interaction mode percentages for the 
remaining background in the nn —>■ search for SKI-IV. 


Systematic Uncertainties 

Systematic uncertainties for the dinucleon decay 
searches are comprised of MC simulation uncertanties 
specific to the signal and background MC, and common 
signal and background uncertainties related to event re¬ 
construction and BDT bias (for modes which use a BDT). 

The largest systematic uncertainties for the dinucleon 
decay signal come from the simulation stage. The domi¬ 
nant systematic uncertainty for all modes come from pion 
final state interaction (FSI) uncertainties. To estimate 
this, we generated separate sets of dinucleon decay MC 
for 24 parameter variations consisting of 1-sigma changes 
in pion-nucleon cross-sections, corresponding to uncer¬ 
tainties determined from fits to pion scattering data. 
This is the same set of variations used in the T2K z/g ap¬ 
pearance result m- The variations affect the cross sec¬ 
tion for both low-energy (charge exchange, absorption, 
quasi-elastic scattering) and high-energy (hadron pro¬ 
duction) interactions. These variations were propagated 
through all stages of the analysis. For the pp —)■ 7r+7r+ 
and pn —>■ modes, each set of MC was processed 

with the trained, unvaried BDT used for the data and 
analysis MC. The systematic error is conservatively taken 
to be the largest percentage change in efficiency across all 
variations upon applying all of the analysis selection cri¬ 
teria. Other significant simulation systematic uncertain¬ 
ties come from Fermi motion (~20%) and correlated de- 
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cay (^10%). These are not passed through the BDT, but 
estimated through a simple re-weighting method of the 
final event sample. For Fermi motion, the re-weighting 
comes from differences between the Fermi momentum 
simulation for dinucleon decay and the Fermi gas model. 
For correlated decay, the effect is not well understood, 
and so a conservative 100% uncertainty is placed on the 
detection efficiency of events with a correlated decay, 
which is 10% at the simulation stage and still about 10% 
in the final event sample. 

The largest overall contributions to background sys¬ 
tematic uncertainties are related to neutrino interactions. 
In particular, for pp —>■ 7r“'"7r+ and pn —>■ the un¬ 

certainties associated with charged current single pion 
production are dominant. This is mainly due to a 40% 
uncertainty in the /tt'^ production ratio. Other simu¬ 
lation systematic uncertainties due to pion FSI and neu¬ 
trino flux are also considered. 




pp - 



Signal (%) 

SK-I 

SK-II 

SK-III 

SK-IV 

Simulation 

35.2 

35.1 

33.6 

38.5 

Reconstruction 

6.0 

8.6 

4.0 

3.2 

BDT 

3.6 

2.2 

4.4 

2.0 

Total 

35.9 

36.2 

34.1 

38.7 

Background (%) 

SK-I 

SK-II 

SK-III 

SK-IV 

Simulation 

29.1 

29.1 

35.8 

26.5 

Reconstruction 

6.1 

8.1 

4.1 

3.2 

BDT 

6.8 

1.0 

4.3 

1.4 

Total 

30.5 

30.3 

36.4 

26.8 


TABLE XVI. Systematic uncertainties for the pp —>■ 
mode. 




nn —^ 

C 

c 


Signal (%) 

SK-I 

SK-II 

SK-III 

SK-IV 

Simulation 

31.1 

34.4 

37.3 

33.1 

Reconstruction 

1.5 

1.7 

4.0 

3.6 

Total 

31.2 

34.4 

37.6 

33.3 

Background (%) 

SK-I 

SK-II 

SK-III 

SK-IV 

Simulation 

13.6 

15.5 

14.5 

13.9 

Reconstruction 

10.9 

18.1 

28.9 

24.3 

Total 

17.5 

24.0 

32.3 

28.0 


TABLE XVIII. Systematic uncertainties for the nn —>■ 
mode. 


1.7%, 2.7%, and 2.3% for SK-I, II, III, and IV, respec¬ 
tively [5] |30]. This uncertainty is propagated into all 
energy-dependent variables, and is propagated through 
the trained, unvaried BDT used for the data and analysis 
MC for pp —>■ 7r+7r+ and pn —>■ The energy scale 

uncertainty is the dominant systematic uncertainty for 
the atmospheric neutrino background in the nn —> 
search. Other reconstruction-related uncertainties taken 
into account are ring separation, particle identification, 
and the ring directional resolution. 

Finally, for the multivariate method, an additional sys¬ 
tematic uncertainty for the training bias of the BDT is 
considered. This is estimated by assessing the difference 
in the testing and analysis distributions of the BDT out¬ 
put for both signal and background. Tables pCVI|XVin| 
break down the systematic uncertainties for each mode 
and each SK period. 




pn ^ 

7r’'"7r‘^ 


Signal (%) 

SK-I 

SK-II 

SK-III 

SK-IV 

Simulation 

33.3 

32.2 

28.4 

35.0 

Reconstruction 

3.3 

1.7 

5.6 

5.6 

BDT 

<1 

1.6 

<1 

<1 

Total 

33.4 

32.3 

28.9 

35.3 

Background (%) 

SK-I 

SK-II 

SK-III 

SK-IV 

Simulation 

22.1 

19.9 

24.0 

27.8 

Reconstruction 

1.8 

1.8 

3.3 

3.8 

BDT 

6.3 

7.4 

10.3 

11.3 

Total 

23.1 

21.3 

26.3 

28.6 


VI. RESULTS AND DISCUSSION 


Since no significant excess was observed in the data, 
we set a lower limit on the lifetime of each dinucleon de¬ 
cay mode. This is done using a Bayesian method, which 
incorporate systematic uncertainties as follows [5^ : 


P(F|n) = A 


e-(rAe+b ) 

n! 


TABLE XVII. Systematic uncertainties for the pn —>■ tt+tt® 
mode. 

The most important reconstruction systematic uncer¬ 
tainty for both dinucleon decay and atmospheric neutri¬ 
nos comes from the energy scale. The energy scale un¬ 
certainty is determined from comparing the absolute mo¬ 
menta of data and MC for several control samples, such 
as cosmic ray muons, their resulting decay electrons, and 
TT^’s from neutral current interactions. Time variation 
in the reconstructed momenta of the control samples is 
also taken into account in the uncertainty calculation. 
The energy scale uncertainty is estimated to be 1.1%, 


xP (F) P (A) P (e) P (b) dXdedb, 

where F is the true dinucleon decay rate, n is the number 
of candidate events, A is the exposure, e is the efficiency, b 
is the background rate, and A is a normalization constant 
that ensures P (F|n) dF = 1. 

The distributions P (A), P (e), and P (6) are taken to 
be Gaussian, while P(F) is uniform for F > 0. The un¬ 
certainty on A is assumed to be negligible. The exposure 
is ^®0-years, not nucleon-years as in most proton decay 
analyses. This is consistent with the exposure used in 
[TT] and [12], which used total and ®®Fe nuclei, re¬ 
spectively. 
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mode 

Frejus limit (®®Fe) This analysis (^®0) 

pp —>■ 7r’'"7r^ 

7.0 X 10^® yrs 

7.22 X 10®® yrs 

pn —>■ 

2.0 X 10®° yrs 

1.70 X 10®® yrs 

nn —>■ 7r®7r° 

3.4 X 10®° yrs 

4.04 X 10®® yrs 


TABLE XIX. Limits on each mode, compared with previous 
results in m- 

Finally, the confidence level is calculated as: 

limit 

C.L.= / p(r|n)dr 
Jo 

and the lower lifetime limit is just the reciprocal of 
Table IXIXl shows the lifetime limits at 90% confidence 
level, compared with the previous limits set by the Frejus 
collaboration in ref. HU. 

The interpretation of our results depends on the un¬ 
derlying physics model. For instance, our limit for nn —> 
TT^TT^ is the most stringent experimental limit for dinu¬ 
cleon decay to light mesons, including m, which mea¬ 
sured Tpp^x+K+ > 1-77 X 10^^ years. However, [11] was 
motivated by supersymmetric models that forbid cou¬ 
plings between quarks of the same generation. If, as with 
the models considered in |33] , such couplings are not for¬ 
bidden, and no preference is given for couplings to first- 
or second-generation quarks, our result for nn 7r°7r° 
provides the strongest general constraint on dinucleon 
decay to mesons. 

As stated in Section I, there is a connection between 
neutron-antineutron oscillation (n-nbar) and dinucleon 
decay to pions, if the same operator is dominant for both 
processes. This is discussed further in Appendix A. Our 
results provide a useful cross-check to n-nbar searches. 
By isolating a single final state for each mode, as opposed 
to the many non-negligible final states of n-nbar [55] . 
these analyses are simple and robust. 


ACKNOWLEDGEMENTS 

We gratefully acknowledge the cooperation of the 
Kamioka Mining and Smelting Company. The Super- 
Kamiokande experiment has been built and operated 
from funding by the Japanese Ministry of Education, 
Culture, Sports, Science and Technology, the United 
States Department of Energy, and the U.S. National Sci¬ 
ence Foundation. This work was partially supported by 
the Research Foundation of Korea (BK21 and KNRC), 
the Korean Ministry of Science and Technology, the Na¬ 
tional Science Foundation of China, the European Union 
EP7 (DS laguna-lbno PN-284518 and ITN invisibles GA- 
2011-289442), the National Science and Engineering Re¬ 
search Council (NSERC) of Canada, and the Scinet and 
West-grid consortia of Compute Canada. 


Appendix: A: Dinucleon decay lifetimes 

Here we discuss the lifetime formula for dinucleon de¬ 
cay to pions. We assume that the same operator which 
mediates neutron-antineutron oscillation (n-nbar) also 
mediates dinucleon decay to pions, which follows from 
crossing symmetry. 

Reference |33| discusses both n-nbar and dinucleon de¬ 
cay with the same operator. This operator takes the 
general form: 


Oab=2 - 

^^sup 

where Oqq is a six quark operator containing two up-type 
and four down-type quarks, and Asup is a suppression 
scale, with the fifth power following from dimensional 
analysis. Indices related to quark generation, color, etc., 
are suppressed. For all first-generation quarks, this op¬ 
erator leads to both n-nbar and dinucleon decay. The 
(free) n-nbar lifetime is given by: 


VII. CONCLUSION 

We have performed searches for dinucleon decay to pi¬ 
ons via the processes pp — >■ ^ pn — >■ and 

nn — >■ 7r°7r° with 282.1 kiloton-years of data with the 
Super-Kamiokande detector. This is the first search for 
dinucleon decay to pions done in a water Cherenkov de¬ 
tector. We do not observe any significant event excess 
for any of the modes, and the data remain consistent 
with the atmospheric neutrino background. Therefore, 
partial lifetime limits of Tpp_,,^+,n.+ > 7.22 x 10^^ years, 
Tp„_,,^+^o > 1.70x10^^ years, and > 4.04x10^^ 

years per nucleus have been set. The limits per nu¬ 
cleus are about two orders of magnitude more stringent 
than previous limits set in m, and constrain models that 
predict AB = 2 processes. 


A^ 

sup 

Te ■ 
^^QCD 


(A.l) 


where Aqc d is a hadronic matrix element near the QCD 
scale, estimated to be 200 MeV. The dinucleon decay 
lifetime tmn is given by [551134] 


S’Pn 


(A.2) 


where is the nucleon mass and pN ~ 0.25 fm is 


the average nuclear density. From equations A.l and A.2 


tnn 


8pN 


A 


2 2 
QC D^nn 


= n 


(A.3) 
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Equation A.3 is mathematically similar to the formal¬ 
ism for comparing free and bound n-nbar lifetimes, al¬ 
though in the latter case the factor T^uc comes from dif¬ 
ferent physics related to the difference in neutron and 
antineutron nuclear potentials. For dinucleon decay, 
Tnuc ~ 1-1 X 10^® while for n-nbar T„„c is two 

orders of magnitude less, with [3^ calculating it to be 


0.517 X 10^^ The current lower limit on the n- 

nbar lifetime, from Super-Kamiokande, is Tnn > 2.7 x 10® 
s [53] • This translates to a lower lifetime for dinucleon 
decay to pions of ~ 2.5 x lO®'^ years, about two 

orders of magnitude beyond the limits obtained in this 
paper. Thus, our results do not place any new constraints 
on n-nbar. 
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